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This article provides a general study of the Hamiltonian stability and the hyperbolicity of vector 
field models involving both a general function of the Faraday tensor and its dual, f(F 2 ,FF), as 
well as a Proca potential for the vector field, V(A 2 ). In particular it is demonstrated that theories 
involving only f(F 2 ) do not satisfy the hyperbolicity conditions. It is then shown that in this class 
of models, the cosmological dynamics always dilutes the vector field. In the case of a nonminimal 
coupling to gravity, it is established that theories involving Rf(A 2 ) or Rf(F 2 ) are generically patho- 
logic. To finish, we exhibit a model where the vector field is not diluted during the cosmological 
evolution, because of a nonminimal vector field-curvature coupling which maintains second-order 
field equations. The relevance of such models for cosmology is discussed. 
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I. INTRODUCTION 

Inflation [1] is usually invoked to explain the isotropy 
and homogeneity of our universe. In particular it has 
been demonstrated that if the dynamics of the universe 
during inflation is dominated by a scalar field, any pri- 
mordial spatial anisotropy is washed out, both at the 
background level Q and perturbation level 0, 0, [H • Sev- 
eral features of the cosmic microwave background (CMB) 
temperature anisotropics seem however not to be fully 
consistent with this prediction. This includes @ the 
low quadrupole (although its statistical relevance is ques- 
tionable), the alignment of the lowest multipoles and an 
asymmetry in power between the northern and southern 
hemispheres. 

It has been suggested that this may be related to 
an early anisotropic expansion during the inflationary 
phase Q- In such a case, it can only lead to an ob- 
servable anisotropy in the CMB at the largest angular 
scales at the price of a fine tuning on the number of e- 
folds during inflation [1, 0, IE]- A natural extension of 
such an anisotropic expansion is to introduce other mat- 
ter fields, besides the inflaton, having the property to 
source the shear. This is the case of vector fields || I^JlfJ] , 
2- forms UM or axions 1121. 



However, vector fields are usually diluted by the cos- 
mological expansion, both during inflation and the mat- 
ter era. Indeed, in a Friedmann-Lemaitre spacetime, with 
metric 1 



ds 2 



-dt 2 



a 2 (t)^iijdx % dx\ 



(1) 



t being the cosmic time, a the scale factor and jij the 
comoving spatial metric, the spatial homogeneity implies 
that the only nonvanishing component of the Faraday 
tensor is Foi- The Maxwell equation reduces to At + 
HAi — [see Sec. Ill B "31 below for a detailed discussion]. 
Thus A 2 cx t 2 ~ 4p if the scale factor scales as t p and A 2 — > 
during the matter era (j> = 2/3) and during inflation 
(p>l). 

This well known fact led to the conclusion that in or- 
der to construct inflationary models driven by a vector 
field, and even to have a slow-rolling vector field during 
inflation, one needs to include either a potential to the 
vector field [3, EH, EH] or a nonminimal coupling [I?], EH . 
The stab ility of these models is actually an ongoing de- 
bate [lM Ha HH IH, HI] . Most of these models have been 
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Vector fields are thus central ingredients in various cos- 
mological models for both the inflationary era and the re- 
cent acceleration. Needless to recall that they also play 
a key role in various extensions of general relativity, with 
the vector-tensor theories [HI, I34l.l35i and more recently 
the tensor-vector-scalar theory [36j that aims at repro- 
ducing the MOND phenomenology, although they have 
several theoretical and experimental difficulties [§3, . 

The goal of this article is twofold. First, we want to 
revisit the dynamics of vector fields during inflation and 
take the opportunity to clarify the structure of theories 
with nonminimally coupled vector fields. A fundamen- 
tal theory should satisfy two necessary conditions: the 
boundedness by below of its Hamiltonian 2 (otherwise the 
theory is unstable [1^]), and the hyperbolicity of the field 
equations (so that the Cauchy problem is well posed fiol]). 
We will derive below the implications of these two con- 
ditions on the vector-field theories we will consider. Of 
course, as soon as these theories are assumed to be ef- 
fective ones, then such conditions need to be satisfied 
only in their domain of validity, but this is still quite 
constraining. 

Section [IT] starts by analyzing theories with a mini- 
mally coupled vector field and a quadratic kinetic term, 
allowing for a Proca potential, and focuses in a second 
part on nonlinear functions f(F 2 ,FF) of the Faraday 
tensor and its dual. We then consider different classes of 
nonminimally coupled theories in Sec. IIIIl To finish, we 
emphasize in Sec. IIVI that there still exist models which 
allow a vector field to be slow-rolling, hence offering an 
interesting cosmological phenomenology. 

Before we start, let us stress that our analysis restricts 
to cases where the vector field is not of constant norm, 
and we refer to Ref. [4l[ where such a case was inves- 
tigated in depth. Let us also stress that the Hamilto- 
nian analysis is more powerful than a perturbative anal- 
ysis around a particular background since the latter can 
only demonstrate the local stability or instability. Hence 
our analysis will generalize in many ways some recent 
results [H M, \m, E3, H concerning the stability of 
vector-field models. 



II. MINIMALLY COUPLED THEORIES 

A. Lagrangian and equations of motion 

As a starting point, let us consider a minimally coupled 
vector field, whose kinetic term is quadratic in its first 
derivatives, and including a potential V(A 2 ), where A 2 = 



A^A^. The most general kinetic term a priori includes 
a linear combination of (V^XVM"), (V^XVA 1 ), 
and (V^A^) 2 . However, the last term can be integrated 
by parts as 

yVaV^V^) 2 = J ^x^—g^M^A^) 
+R^A„A 1/ ], (2) 

so that only a linear combination of the first two terms 
needs to be considered in flat spacetime. However, in 
curved spacetime, the extra term A^A V is a particu- 
lar nonminimal coupling to gravity. 

Let us first recall that, in flat spacetime, the only 
ghost-free vector theory in the above class is the stan- 
dard Maxwell Lagrangian (called Proca Lagrangian in 
the massive case [42| ) 

^Maxwell = — 4^ 2 ' 00 

where F 2 = F 2 „, and F^ v = d fl A l/ — d v A^ is the Faraday 
tensor. Indeed, if we consider a Lagrangian 

C = a(dM 2 + /3(a M ^,)(0M") - V(A 2 ), (4) 

we deduce from F^i = At — dtA® that the conjugate mo- 
menta 7r M = dC/dA^ read 

tt° = 2(a + 0)A O , tt* = -2aA % - 2f3d l A a . (5) 

If a + (3 ^ 0, the field A$ is thus dynamical. This can 
also be illustrated by writing the Euler-Lagrange equa- 
tion deriving from (U) 

a\JA" + fid^^ -V'A", (6) 

together with its divergence 

{a+P)U<t> = -d x {V'A x ), (7) 

where = d^A** and V = dV7d(,4 2 ). Although the re- 
placement of the derivative d^A* 1 by a scalar field would 
be illicit 3 in the Lagrangian (|4|), one may do so in the 
field equations, and Eqs. ©-([7]) show that the model 
describes a transverse vector field (A^ with d^A^ = 0) 
together with a scalar degree of freedom (j>. These equa- 
tions also underline that some degrees of freedom become 
nondynamical when either a = or a + [3 = 0, as will be 
discussed below. 

Let us first consider the generic case where a/0 and 
a + /3^0. Then the Hamiltonian density 7i = n^A^ — C 



2 More precisely, the spatial integral of the Hamiltonian density 
over any localized state should be bounded by below. Since 
such localized states may be constructed from a superposition 
of sinusoids, at least at linear order, one may also compute the 
Hamiltonian density for such spatial sinusoids. 



Redefining a derivative as a fundamental field in a Lagrangian 
obviously loses some dynamics, as illustrated by the trivial case 
of a scalar-field kinetic term C = — (c9 M </?) 2 , which would give an 
adynamical vector C = — V 2 if one redefined V M = dptp. 
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takes the form 

(vr ) 2 (tt* + 2/3c^ ) 2 



H = 



a(cU ) 2 



4(q + /3) 4a 



Since 7r l and OiAq are independent from each other, and 
since the quadratic form —x 2 /a + ay 2 is not positive 
definitive (whatever the sign of a), we conclude that Tt 
can take arbitrary large and negative values, and thereby 
that the theory is unstable. This is the well-known re- 
sult that the massive vector contains three modes of 
positive energy, but also an extra helicity-0 ghost. 

On the other hand, if a + /3 — 0, which corresponds 
to the usual Maxwell Lagrangian ([3]), then Eq. ([5]) yields 
the primary constraint ir° = and the scalar mode is 
no longer dynamical. We will recall in Eq. (fT8|) below 
the standard result that — a — (3 > and V' > are 
necessary conditions for the Hamiltonian to be bounded 
by below. 

The other particular case for which expression ([5]) for 
the Hamiltonian cannot be used is when a = 0. After 
integration by parts, this corresponds to a simple kinetic 
term of the form (3(d^ L A^) 2 . The conjugate momenta read 
then tt = 2(3(A - d t A l ) and tt 4 = 0, so that only the 
helicity-0 degree of freedom contained in the vector A^ is 
now dynamical. This case should thus be considered as a 
scalar theory rather than a vector one (although it differs 
from standard scalar theories because of the secondary 
constraint 2AiV — 8^° imposed by the field equations). 
We will thus not consider it any longer in this paper. Let 
us just mention that the first term contributing to the 
Hamiltonian density Tt = (tt° ' + 2[3d l A 1 ) 2 / '4/3- (3{d l A 1 ) 2 + 
V can obviously be made positive by choosing /3 > 0, 
but that this does not suffice to guarantee the stability 
of the model because the second term is then negative. 
The fact that it is not independent from 7r° complicates 
the analysis, but when V' — const. > 0, for instance, 
it is easy to build consistent initial conditions such that 
Tt — > — oo, thereby proving that the model is unstable in 
such a case. 

This analysis underlines that vector-field theories are 
generically unstable when their kinetic term does not re- 
spect the gauge invariance A^ — ► A^ + d^X, because the 
Aq component is then a ghost degree of freedom 4 . Before 
studying nonlinear vector actions, let us underline that 
the above Hamiltonian analysis is fully changed in the 



4 Let us recall that a ghost is defined as a field with negative kinetic 
energy, not to be confused with a tachyon, a field with negative 
mass squared, or more generally with a potential which is un- 
bounded by below. Both cases correspond to unstable models, 
but tachyons involve a time scale whereas the presence of ghosts 
implies an instantaneous disintegration of vacuum in quantum 
mechanics |39(. Note that a field may be both a tachyon and 
a ghost, but that the corresponding model is then even more 
unstable. 



case of a constant-norm vector field; see Ref. 41] for a 
detailed analysis of this interesting case. 



B. Function of F 

Let us thus consider now nonlinear functions of F 2 , 
i.e., gauge-invariant kinetic terms by construction, in La- 
grangians of the form 



£=-f(F 2 )-V(A 2 ). 



(9) 



The associated field equation for the vector field is then 
simply given by 



1. 



V M (f'Fn = -V'A 



(10) 



where a prime denotes a derivative with respect to the 
argument of the function, namely /' = d//d(F 2 ) and as 
before V = dV/d(A 2 ). Note that /' should never vanish 
otherwise the Cauchy problem would be ill-posed. From 
the definition of the Faraday tensor, we always have 

d a F^ + d^F va + 8 v F a ^ = 0, (11) 

and the divergence of Eq. (JTUJ) implies 

V V {V'A V ) = 0. (12) 

When V' ^ 0, this is an extra constraint that arises from 
the fact that the action is no more invariant under A^ — » 
Afj, + <9 M A, even if the kinetic term independently is. 

1. Hamiltonian analysis 

Since we have F 2 = F-j—2F§ i in Minkowski spacetime, 
the conjugate momenta read 

7T° = 0, (13) 

which is a primary constraint, and 

vr* = 4/'(F 2 ) x (Ai - d.Ao) = 4/'F 0l . (14) 

Since Aq does not appear in Lagrangian @, Aq is an 
auxiliary field. This means that the field equation for 
Aq involves no time derivatives and can be used as a 
constraint that eliminates a field variable, in the case at 
hand Aq. 

The Hamiltonian density is thus given by 

H=^.+Tr*d l Ao + f(F 2 ) + V(A 2 ). (15) 

The Aq dependency of Tt can be eliminated by first per- 
forming an integration by parts (in which ^diAg be- 
comes — Aodiir 1 ) and then using the secondary constraint 
[7To,W] = 0. This secondary constraint ensures that the 
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primary constraint (|13[) is consistent with the equations 
of motion, and it takes the form 5 

cW = -2V'A . (16) 

Actually, it turns out to be the Euler-Lagrange equation 
(IIOP for v = 0, rewritten in terms of conjugate momenta, 
and it reduces to the Gauss law when V' = 0. Note 
that, in general, there may be further constraints arising 
from the consistency of the secondary constraints with 
the equation of motion, and so on. The distinction be- 
tween primary and secondary is not important and they 
are just constraints that we consider on the same footing. 
It follows that 

n = ^ J + { -^ + f(F 2 ) + V(A 2 ), (17) 

if we assume that V' 7^ 0. [In the case where V' = 0, 
then diW 1 = from Eq. (fTS)) . so that Tt does not in- 
volve any term oc (diir 1 ) 2 .] This is a function of the field 
Ai, its spatial derivatives diAj, its conjugate momentum 
7T 4 and its derivatives diir J , since the argument of the 
function / can be expressed as F 2 = F 2 j — 27r 2 /(4/') 2 
and Aq can be eliminated by resolving Eq. (|16[) ; hence 

Equation (jTTJ) shows that it is necessary that /' be 
positive for Ti to be bounded by below. Indeed, if there 
existed a value, say F 2 , where f'(F 2 ) < 0, then one could 
construct initial conditions where irf — > 00 and i 7 ! 2 - — * 00 
while keeping F 2 = F?. - 2?r 2 /(4/') 2 constant. The first 
term of the r.h.s. of Eq. (fT7| would then tend towards 
— 00 whereas the other ones would remain finite. 

Similarly, V must also be positive for Ti. to be bounded 
by below. Indeed, using the secondary constraint (Ti7i)) . 
the contribution of the potential to Eq. (fT?"]) reads 
(d l n 1 ) 2 /2V + V = 2A 2 V'(_A 2 ) + V{A 2 ). If there existed a 
value, say A 2 , where V'(A 2 ) < 0, then one could choose 
initial conditions where Aq — > 00 and Af — > 00 while 
keeping A 2 = A 2 — Aq constant, and the Hamiltonian 
would thus diverge towards —00. 

On the other hand, note that the potential V itself 
does not need to be bounded by below, contrary to what 
one may naively believe from Eq. (| 1 T|) . Indeed, the pos- 
itive contribution 2AqV' (A 2 ) can compensate negative 
ones coming from V(A 2 ). For instance, for a mono- 
mial V(A 2 ) = k(A 2 ) n , where k and n are constants, the 



5 Note that the constraint {16} will be general for any theory in 
which 7r° = and 

dC dC 
ddiAo ~ ~aV 

since then the Euler-Lagrange equation implies 



This is the case of all theories in which the kinetic term of the 
vector field involves only functions of F 2 and FF (see Sec. Ill CI 
below) . 



contribution of the potential to the Hamiltonian reads 
k [(2n - I) A 2 + A 2 ] (A 2 ) n -\ therefore it is bounded by 
below if k > and n is a positive odd integer. In 
such a case, V — fcn( J 4 2 )™ -1 is consistently positive, 
but not V itself since it can have any sign. The par- 
ticular case n = 1 corresponds to the standard massive 
Proca field, with V = \m 2 A 2 , i.e., IV 1 = m 2 > 0. Then 

V = — ^(diW 1 ) 2 /m 2 + \rr? Af contains a negative term 
which can blow up for some specific initial conditions, 
but it is counterbalanced by the second term of (|17jl . 
+ (diir 1 ) 2 /m 2 . The above example of a monomial also 
illustrates that V' > is not a sufficient condition. In- 
deed, if one chose k < and n odd and negative, then 

V would always be positive but TL would diverge to- 
wards —00 for initial conditions such that diir 1 = and 
A 2 ^00. 

Some negative contributions coming from f(F 2 ) may 
also be compensated by irf/Af. This is again what hap- 
pens in the massive Proca (or pure electromagnetic) case, 
where f(F 2 ) = F 2 /4 = F 2 - tt 2 /2 but tt 2 /4/' = tt 2 , so 
that 

is clearly positive. 

Since there is no obvious necessary and sufficient con- 
ditions warranting that the Hamiltonian (|17[) is bounded 
by below in the most general case, this should be checked 
explicitly for any specific theory at hand, recalling that 
compensations between terms often occur. 

2. Hyperbolicity 

The second necessary condition that a field theory ^ 
should satisfy, is that its field equations (JTOJ) are hy- 
perbolic, i.e., that their second derivatives are of the 
form G^d^dv, with G^ v an effective metric of signa- 
ture — h ++ (its timclike direction, corresponding to the 
negative eigenvalue, should also be consistent with the 
standard time direction of g^)- These second deriva- 
tives can be written as an operator acting on the vector 
field A MJ 

[/' x - d a d v ) + Af'F^FP^dp] A°. (19) 

Our first difficulty, with respect to the better studied 
case of scalar "k-essence" Lagrangians [13, EE 13 is 
that A^ has four components and that the above op- 
erator is not diagonal. In order to diagonalize it, it is 
convenient to first remove the —f'dad 1 ' contribution to 
Eq. (|19p by fixing the Lorenz gauge, namely by adding 
X(d fl A t1 ) 2 to Lagrangian ([9]), where A is a Lagrange mul- 
tiplier. In this gauge, the operator (TT9")) becomes of the 
form /'□! + 4f"\v){v\, where the Dirac ket \v) repre- 
sents F^ v d v . One finds thus immediately that its four 
eigenvalues (still as an operator) are three times /'□, and 
once f'n+4f"{v\v) = fD + ifF^F^d,,, this fourth 
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"eigenoperator" acting in the direction of \v). Obviously, 
the operator /'□ is hyperbolic of signature — h++ if and 
only if 

/' > 0. (20) 

The fourth eigenoperator may be written as G pl, d^d v , 
where G pv = f g pv + Af"F^F v p is an effective metric 
in which the fourth component of the vector (in our 
specific diagonalizing basis) propagates. The simplest 
way to analyze its hyperbolicity, and to ensure that its 
timelike direction is consistent with the one of g pu , is to 
diagonalize the matrix G pp g pv and impose that its four 
eigenvalues are positive. [Note that our analysis uses two 
different diagonalizations: first a 4 x 4 matrix, with oper- 
ator values, acting on the vector A^; now the quadratic 
differential operator G^d^dv, acting on one particular 
component of A^. ft happens that G pp g pu is again a 
4x4 diagonalizable matrix.] These eigenvalues read 

/' + f"t% ± f"^(F*,) 2 +(F tu ,F>»')\ (21) 
where F^ v is the dual of the Faraday tensor, 

Ffiu — 2^ iivpoF 9 7 (22) 

Sfifpa being the totally antisymmetric Levi-Civita tensor 
such that £0123 = +1. An elegant way to derive these 
eigenvalues is to separate F^ into standard electric (E p ) 
and magnetic (B p ) field contributions according to an 
observer with unit velocity w M . Then, in the generic case 
where E and B are not parallel, one may study the action 
of the operator F p p F p a on the four linearly independent 
vectors E p , B p , u p and g^e a ^ s u a E^ B 5 , and one finds 
that the spaces spanned by the first two and the last two 
are stable under the action of this operator. In other 
words, its matrix is constituted of two 2x2 blocks. Its 
eigenvalues are then easy to compute, and they happen to 
be the same for each block. [The particular cases where 
E and B are parallel or one of them vanishes are easier 
to study along the same lines, and one can check that the 
result (|21[) remains valid.] 

The simultaneous positivity of eigenvalues (f2Tj) im- 
poses thus the subtle second condition 

FlJ" + f > \f"\\j{F^f+(F i ,uF^)\ (23) 

When F 11V F IJ ' V = 0, i.e., when the electric and magnetic 
fields are orthogonal, this inequality imposes both /' > 
[already necessary in Eq. (f2"0"|) above] and 2F 2 f" + /' > 
0. This should be compared to the case of scalar k- 
essence models, whose Lagrangians are functions f(s) of 
the standard kinetic term s = (d^ip) 2 . Then the hyper- 
bolicity of the field equations implies both /' > and 
2s/" + /'> 0,111 Si. 

The fact that the inequality (12"3")) depends on two inde- 
pendent relativistic invariants constructed from the elec- 
tric and magnetic fields, namely F? v — 2(B 2 — E 2 ) 



and F pv F pv = — 4E • B, underlines that it should al- 
ways be possible to violate it by choosing appropriate 
initial conditions on a Cauchy surface. For instance, if 
/"(0) ^ 0, then one may choose a configuration where 
F i = — Fio = F 2 s = —F 32 and all other components 
vanish. Then F 2 vanishes whereas E • B can be chosen 
as large as one wishes. This suffices to violate inequal- 
ity (f2"3")l . and thereby to prove that the field equations 
cannot remain hyperbolic in all physical situations. If 
the considered theory is such that /"(0) = 0, i.e., that 
its Lagrangian does not contain any term proportional 
to (F 2 V ) 2 , then we need to refine slightly the reasoning: 
We choose a value of F 2 such that f{F 2 ) ^ and we 
add to it a contribution E • B increasing the value of the 
square root in (f2"3"|) . while keeping F 2 = 2B 2 — E 2 con- 
stant. The only possibility to always satisfy inequality 
(p3| would be to assume that f"(F 2 ) = for any F 2 , so 
that f(F 2 ) = kF 2 + 2 A (where k and A are constants) 
would merely describe standard Maxwell (or Proca) the- 
ory plus a cosmological constant. 

In conclusion, although theories (f53|) can have a Hamil- 
tonian (fT7|) bounded by below for specific functions 
f(F 2 ), there always exist situations in which the field 
equations are not hyperbolic, because inequality (f2"3"l) is 
violated. The only safe case is the standard Maxwell La- 
grangian (with an optional Proca potential). Of course, 
if such models are considered as effective theories, then 
all the above conditions must be satisfied only in their 
domain of validity. But if one uses such an effective the- 
ory in situations where Eq. (|23p may be violated, then it 
just loses any meaning, since the Cauchy problem is no 
longer well-posed. 

3. Cosmological dynamics 

Let us investigate the cosmology of the models de- 
scribed by Lagrangian The stress-energy tensor of 
such a vector field is given by 

Tpu = 4/'F A M F A!/ + 2V'A P A V - (/ + V)g^. (24) 

Different roads can then be followed. In particular, it 
is clear that the vector field induces the existence of a 
particular spatial direction, in contradiction with the hy- 
pothesis of isotropy underneath the form {T]) of the met- 
ric. One should then consider anisotropic cosmological 
spacetimes, such as Bianchi universes, which character- 
ize the anisotropy, or try to recover isotropy by invoking 
the existence of N vector fields with random directions 
and similar initial magnitude [l?] ]. 

For the sake of simplicity, we investigate the dynamics 
of a test vector field, the dynamics of which is described 
by Lagrangian ([9]) in a cosmological spacetime with met- 
ric ([I]). We can then always decompose the vector field 
as 

Ap = (A , aBi), A» = (-A , -B l ) , (25) 
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with B l = -f^Bj. In Cartesian coordinates, homogene- 
ity implies that diA^ — so that the only nonvanishing 
component of the Faraday tensor is 



F 0l = At = a(Bi + HBi 



(26) 



where H = a/a denotes the Hubble function. As ex- 
pected, Aq will not enter the equation of evolution and, 
as long as V ^ 0, the field equation (jTOj) implies in Carte- 
sian coordinates that A = and 

/' (F w + HF oi ) + fd (F 2 )F m = —V'At. (27) 

Since F 2 = 2F 0i F 0i = -2F 0j F 0k ^ k /a 2 = -2C t C l = 
— 2C 2 , this equation rewrites as an equation for Bi as 



B, 



3H-2^d t (C 2 ) 



Bi 



{2H 2 + H) 



V f" 

— -2H J -d t (C 2 ) 



Bi = 0, (28) 



where we use that /' should not vanish, or equivalently 
as the system 



Ci + 2 
B 4 + HB 



H-^d t (C 2 ) 



1 V 

Cj = Bj 

2 f 



Ci. 



(29) 
(30) 



In that particular case, we deduce that the energy density 
of the field, pa = —T®, is 

PA = if'C 2 +f + V. (31) 

Note that the isotropic pressure Pa =T-/3 is given by 

Pa = -tf'Cf + %V'B\ -f-V. (32) 



For such a vector field the pressure is however not 
isotropic and there is a contribution of the vector field 
to the anisotropic stress (i.e. the transverse and traceless 
part of the stress-energy tensor) 

ttJ = -4f {c,C l - \c 2 5]^ + 2V \&Bi - \B 2 5pj . 

(33) 

From the expression of the energy density and anisotropic 
stress, we see that, in order for the vector field to play 
any significant role, one needs either Q or Bi not to be 
diluted during the expansion 6 . 



In the standard case of the Maxwell theory (/' = 1/4 
and V — 0), it is obvious that Eq. (f2"9"| implies that 
d oc a~ 2 . We then conclude that pa oc a -4 and the 
vector field energy density is diluted with respect to the 
matter fields driving the expansion of the universe. In- 
deed, this could have been deduced from Eqs. (|3i p - (|3"2"l) 
which imply that, as expected, the equation of state of 
the homogenous fluid is 1/3. 

Again, in the Proca case (/' = 1/4 and V' ^ 0), the 
vector field can play a role if it is not diluted, i.e., if Bi ~ 
const, is a solution of Eq. (|28p. This happens if the 
coefficient of Bi is small compared to H 2 , and the energy 
density pa of the vector field is then almost constant. 
However this requires that V < 0, as initially proposed in 
Ref. [HI , in contradiction with the Hamiltonian analysis 
above. 

In the general case, assuming that inflation is described 
by a de Sitter phase, i.e., H is constant, the solution Bi 
constant can only been reached under the condition that 
2H 2 + V 12 f -2H(f" I f)dt(Cf) < H 2 . This is actually 
impossible since Ci — HBi is also constant and V / f is 
positive. This can be generalized to the case of slow-roll 
inflation for which H = —eH 2 . A configuration with Bi 
constant can be reached if 



2H 2 



XL 

2f> 



+ AeH B 



■r 
r 



0. 



(34) 



Since V / f > 0, this is possible only if /"//' is of order 
1/e and e < 0. Such a fine tuning is very unnatural 
since / enters the vector field sector, while e is set by 
the matter driving the inflationary era. On the other 
hand, a configuration with d constant requires, from 
Eq. dnnj), that aB { = G t f adt. But Eq. implies that 
2Ha = —(V'/2f) J adt, which is impossible as long as 
the universe is expanding. In conclusion the vector field 
cannot play a cosmologically relevant role. 

This is confirmed by a more general argument. Let us 
introduce 



= C 



4' = B 2 



dB l = ^cjjtPp 



(35) 



with p 2 < 1 and 4> and ip positive. From the system (|2"9"|) - 
(l30|) . we can extract the following set of equations de- 



In particular, if one relaxes the hypothesis of isotropy and de- 
scribes the universe by a Bianchi I space-time with metric 

ds 2 = -dt 2 + a 2 (t)7ij(i)dxMx J ', 

the shear Oij = is sourced by this anisotropic stress and 

evolves as <j\ +3HrTj = 8TrGiTj so that a nonvanishing anisotropic 



stress can source the shear which decays as a~ otherwise; see 
Ref. Q. Such a vector field, even if it does not influence the 
dynamics of inflation can be at the origin of an homogeneous 
shear, along the line of Ref. [jj. Note that the evolution of the 
vector field is modified so that equation l|28j l has now a r.h.s. 
2tr^ Cj, and dt(C 2 ) now contains a shear-dependent contribution 
-2<T« CiCj. 
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scribing the relative evolution of Bi and C, : 

V 



f" 



-m<j>- n—y/W), 




In this system, V' is a function of —-0 and /" and /' 
are functions of —2<j) so that the system has been written 
as a dynamical system. Its fixed point, characterized 
by <fi — ip = 0, must be such that Vip/f = — 4<£, in 
contradiction with V'/f > unless <f> = and V'^f = 0. 
Even if V = 0, setting ip = and <f> = in Eq. (I57| 
implies ^ = as soon as _ff 7^ 0. Therefore the unique 
fixed point of this dynamical system corresponds to (f> = 
ip = 0, i.e., to a strictly vanishing vector field. 

In conclusion, slow-rolling solutions can be constructed 
at best via an unnatural fine-tuning, and moreover, these 
solutions are not fixed points of the dynamics. We con- 
clude that such vector fields will be diluted and play no 
role in cosmology. 



C. Introducing FF 

Since the contraction FF = F^F^ appeared in the 
previous hyperbolicity analysis, we are naturally led to 
consider an extension of theory ^ of the form 



£ = -/ \F ,FF\ - V(A Z ). 
In the following, we will set 

X = F 2 and Y = FF, 



(39) 



(40) 



and denote as fx and fy the partial derivatives of F 
with respect to X and Y, respectively. The field equation 
deriving from l|39p can thus be written as 



7 , (/. 



xF^ + fyF^ 



1 



V'A". 



(41) 



Note that in the particular case in which / = FF, this 
equation is empty because of the identity d^F^ = (i.e., 
d 2 ^4 = in Cartan's exterior-derivative notation, namely 
Maxwell's first set of equations, F^ y . p ^ = 0). 

An example of such theories, though it is an effective 
one, is the Euler-Heisenberg corrections [42|, HU to the 
Maxwell Lagrangian ([3]), which take into account the vac- 
uum polarization. It is given by the Lagrangian 



£eh = 



90m! 



(F^F^ 
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(42) 



where a is the fine-structure constant and m e the mass 
of the electron. It is derived formally as the first term 



of an expansion when a 2 — > 0, and its domain of valid- 
ity is precisely when such nonlinear corrections remain 
small with respect to the standard Maxwell theory ((3|. 
In this domain of validity, the Hamiltonian density is 
positive and the field equations are hyperbolic, therefore 
none of the following discussions need to be done. On 
the other hand, as soon as a Lagrangian of form (f3"9"l) 
is considered as defining a fundamental theory, or when 
one wishes to study its predictions in a domain where 
nonlinear effects are significant, then both the stability 
and the well-posedness of the Cauchy problem need to 
be analyzed carefully 



1. Hamiltonian analysis 

As in the previous sections, we need to compute the 
Hamiltonian density and we restrict to a Minkowski back- 
ground spacetime. The two relativistic invariants J40|) 
reduce to 



X = F 2 -2F 2 and Y = 2e ijk F 0i F jk , 



(43) 



where we have set e^fe = eoijfc. It follows that the conju- 
gate momenta take the form 



ir Q = 0, 

7T* = 4f x Foi-2f Y e ijk F jk , 
and the Hamiltonian density reads 



(44) 
(45) 



H = yj £rS ijk KiF jk + ^d t A + f + V, (46) 

We are here assuming fx 7^ 0, and will consider the 
particular case of functions of Y alone in Sec. Ill C 3l below. 
The field equation (|41 [) reduces, as expected from the 
comment in footnote [51 to 



d.y = -2V'Aq. 



(47) 



Integrating by part the term n l diAo and then using the 
secondary constraint to eliminate Aq, we end up with a 
Hamiltonian density 



n 



4/x 



{di^y f Y 
Tx 



2V 



^mdjAk + f + V. (48) 



[This expression assumes that V ^ 0. When it van- 
ishes, the second term oc (diir 1 ) 2 merely disappears be- 
cause Eq. (|4T[) implies diir 1 = 0.] 

As for the simpler case of functions f(X) considered 
in Sec. IIIB1 one needs to check that the Hamiltonian 
density |48p is bounded by below for each specific model 
one is considering. 

A necessary condition is that fx be positive. Indeed, 
in terms of Ei = Foi and B l = e % ^ k djA kl the Hamiltonian 
may be rewritten as 



H = 4f x E 



2 (di^l 



f-Yfy + V. 



(49) 
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Now, one can let E — > oo while keeping constant the ar- 
guments X = 2(B 2 - E 2 ) and Y = -4E • B of the func- 
tion / and its derivatives. [This can be performed for 
instance by setting E = J X/2 sinhp, B = ^/X/2coshp, 
cos(E, B) = — Y/(2X sinhpcoshp), and letting the pa- 
rameter p — > oo.] Therefore TL could take arbitrary large 
and negative values if we had fx < 0. 

Specific sufficient conditions may also be written to 
ensure that Ti. is bounded by below. For instance, it 
would obviously suffice that fx > 0, V > and both 
/ — Yfy an d V are bounded by below. However, this is 
far from being necessary, since the positive contribution 
coming from 4fxE 2 can compensate a negative one due 
to / - Yf Y , and that (d l ir 1 ) 2 /2V 1 + V may be bounded 
by below even if one of the terms can diverge towards 
— oo. This is what happens in the standard Proca case 
discussed in Eq. (fTgj) above. 



2. Hyperbolicity 

Following the same lines as in Sec. lII B 2l equation (HT]) 
for the propagation of the scalar field can be rewritten as 
an operator acting on A" , 

fx{8 v a U-d a d») 

+4 [fxxF^F^ + fvyF^F^) d^dp 
+4fxr (F^F^ + F^FP^j D^dp. (50) 

For specific particular cases, it is possible to diagonalize 
its action as independent operators acting on the com- 
ponents of A a , and their hyperbolicity can then be ana- 
lyzed as before by working in the generic basis E^, B^, 
it M , g^ >1 t a fi GV u a B v E a . However, the first diagonalization 
is quite involved and we did not derive the most general 
conditions which must be satisfied. Moreover, the anal- 
ysis of necessary or sufficient conditions on f(X,Y) en- 
suring hyperbolicity is also a difficult task. Therefore, we 
merely conclude that for each specific model, one should 
check both the boundedness by below of the Hamiltonian 
density P5|) - (P5]) and that the matrix of operators ([50]) 
defines hyperbolic equations for all physical components 
of the vector A^. However, we shall see in Sec. Ill C 41 
below that this class of models (|39[) does not answer the 
question we are addressing in the present paper, i.e., that 
the vector field is necessarily diluted by the cosmological 
expansion. 

3. Particular case of f(FF) 

The above Hamiltonian analysis assumed that fx^O, 
therefore it cannot be followed in the special case where 
f(FF) does not depend on F 2 . In such a case, it is 
straightforward to show that the corresponding Hamil- 
tonian density is bounded by below only if V' > and 



/ — Yfy is itself bounded by below [the discussion con- 
cerning the potential V is the same as below Eq. (fT7|) ]. 
However, as in Sec. Ill B 21 above, the analysis of the hy- 
perbolicity of the field equations suffices to exclude these 
models. Indeed, the field equations read 

2<9 M (F^/') =A v V, 

that is to say 

2F^ v d tl f' = A V V'. (51) 

This equation already shows that no propagation of per- 
turbations can be defined through a spacetime hyper- 
surface where the background value of F^ v happens to 
vanish. This suffices to underline that this class of mod- 
els is pathological. One may anyway mimic the analysis 
of Sec. Ill B 21 and diagonalize the differential operator 
acting on A^ in Eq. (|5~lj) . One finds that three out of the 
four components do not propagate because they have a 
strictly vanishing differential operator. The fourth com- 
ponent is differentiated by the operator G^d^d^, where 
Qtiu _ 4 jii pwpv^ plays the r61e of an effective me tric 

in which perturbations propagate. The same reasoning 
as in Sec. Ill Bl2"1 above then shows that the eigenvalues of 
the matrix G Mf 'g pv cannot all be simultaneously positive, 
and therefore that this last differential operator is not 
hyperbolic either. Indeed, one would need to satisfy the 
strict inequality 



[F, v ff" > \f"\^(F^) 2 + (P^F, U ) 2 , (52) 
which is impossible. 

4- Cosmological dynamics 

In the particular case of an homogeneous space-time, 
and as detailed in Sec. lIIB~3l the only nonvanishing com- 
ponents of the Faraday tensor are F$i so that only Fij is 
nonvanishing and thus, it implies that FF = 0. 

As a consequence, the field equation (fH) leads to the 
same equation as for the case of a function of X alone 
that is to Eq. (|2"8"|) with the function f(X) replaced by 
f(X, 0). The cosmological dynamics remains unchanged 
and the conclusions of Sec. Ill B 31 are not affected. 

5. Conclusions and remarks 

Our analysis shows that f(F 2 ) models do not satisfy 
the hyperbolicity conditions (unless /' = const.), and 
that one must then extend them to f(F 2 ,FF). This is 
needed if the model is considered as a fundamental the- 
ory, but also in the domain of validity of an effective 
one. As we shall also see below, an interesting cosmo- 
logical phenomenology can generically be obtained only 
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when the nonlinear corrections become comparable to the 
lowest-order F 2 kinetic term. The hyperbolicity condi- 
tions need thus to be satisfied in such a case, even if the 
model is assumed to be effective. 

Independently of these conditions, we also showed that 
the only fixed point of the cosmological dynamics corre- 
sponds to A p — 0, so that the vector field is diluted 
during the cosmological expansion, and therefore cannot 
play any significant cosmological role. 

We could have imagined more complex terms such 
as F^F v PF p ^ or F^F l, PF pa F av . However, one can 
check that the first combination strictly vanishes while 
the second can be rewritten as a function of F 2 and 
FF, so that our analysis above already considered such 
possibilities. Let us also point out that terms such 
as {d f _ L A v ){d"A p ){d p A^) generically excite the helicity-0 
ghost degree of freedom. 

D. Constant norm vector field 



Following the same approach as in Sec. Ill C 11 we get 

H = AfxE 2 - y/A 2 -vV-n + f- Yf Y + V. (58) 

We conclude that whatever the functions / and V, this 
Hamiltonian density is not bounded by below because 
one can always let X, Y and A 2 constant while letting 
V • 7r go to infinity. 

It should be underlined that the above conclusion 
only applies to kinetic terms of the form ([53")) . As 
shown in Ref. :4l|, more general kinetic terms for a 
constant-norm vector field, of the form Ci(d„A u ) 2 + 
c 2 (0p A" ) 2 + c 3 (0 M A, ) {d v A» ) + c 4 ( A» A v ) 2 , can be con- 
sistent for specific ranges of values of the constant coeffi- 
cients Ci 2,3,4, i.e., define stable and well-posed field theo- 
ries and even pass solar-system and binary-pulsar tests of 
relativistic gravity. The same analysis has not yet been 
generalized to nonlinear functions of such kinetic terms, 
nor to variable coefficients (depending on some field). 



Given the conclusion of the previous analysis it is in- 
teresting to consider similar theories but with the con- 
straint that the vector field has a constant norm. Gen- 
eral study of constant norm vector fields have been dis- 
cussed notably in Ref. [4l[, and they play an important 
role for instance in the construction of MOND-inspired 
theories [H,^. 

We may thus consider Lagrangians of the form 



C 



f*,ff) 



V(A 2 ) + X(A 2 -v), 



(53) 



where A is a Lagrange multiplier and v a number. The 
extremization of the action with respect to A gives the 
constraint 



A 2 



(54) 



and A^ is timelike (resp. spacelike) when v < (resp. 
v > 0). The norm-fixing term does not change the ex- 
pression of the conjugate momenta which are still given 
by Eqs. (|44[) - ([4"S"|) . The equation of motion gets an extra 
term 



d t ir l = -2V'A + 2\A 



(55) 



It cannot be used to eliminate Aq from the Hamiltonian 
density since it is now used to fix the value of A. Instead, 
we use Eq. ([54]) to get 



A = Jj%-v. (56) 
We conclude that the Hamiltonian density simplifies to 

n = ^--Ja 2 - vidi^y^e^mdjAk+f+v, (57) 

4/x v JX 



where the only difference with the expression (|48p lies in 
the second term. This expression should be compared to 
the result of Refs. [13, HJ. 



III. NONMINIMAL COUPLINGS 

The results of the previous section drive us to consider 
theories with a standard kinetic term. This section fo- 
cuses on models satisfying this constraints but involving 
a nonminimal coupling to gravity. This class of models 
is of particular interest in cosmology because it has been 
argued that when such a coupling exists the vector can be 
slow-rolling [17[ and the stability of this models has been 
debated with different conclusions [H, [HI, H HI . 

We already saw, in Eq. @ above, that nonminimal 
vector-metric couplings of the form ' A^A V are gener- 
ated by a mere integration by parts of a general vector 
kinetic term in curved spacetime. Such a term, together 
with a RA 2 coupling, has been considered in chapter 5.4 
of Ref. In the following, we will not study R^A^A,,, 
whose mathematical and phenomenological consequences 
are similar to those of RA 2 . However, we will consider 
the more general case of nonlinear couplings to a func- 
tion of A 2 in Sec. IIII A[ and show that the corresponding 
models are unstable. We will also consider couplings to 
a function of the Faraday tensor F in Sec. IIIIB1 but un- 
derline that instabilities are also generic in such a case. 



A. A case 



Jordan frame 



Let us first consider models of the class 



S = J A'xJ—g 



R ^(A 2 )--F 2 -V{A 2 ) 



2k v ' 4 

+'S'matter [V'm I 9 1 



(59) 



where k = 8irG, g M „ denotes the Jordan frame metric, 
and we define F 2 = F^F^g^g 1 "' and A 2 = A^A^ . 
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is an arbitrary positive function and the particular case 
"J = 1 + 8irG£,A 2 has been extensively studied in the lit- 
erature G is the bare gravitational constant. It is 
not the constant that would be measured in a Cavendish 
experiment since the vector field is responsible for an in- 
teraction. As in the case of scalar-tensor theories [It}, the 
Jordan metric is the "physical metric" since the matter 
fields are universally coupled to g^ . This metric defines 
the lengths and times actually measured by laboratory 
rods and clocks, since they are made of matter. All ex- 
perimental data have their usual interpretation in this 
frame. 

The equation of motions, obtained by variation with 
respect to the vector field, is given by 



V„F'" - ( 2V--V' ) A" = (i 



(60) 



which generalizes the Maxwell equation. As previously, 
a prime denotes a derivative with respect to the argu- 
ment, V = dV(X)/dX . The divergence of this equation 
implies that 



2V' 



R 



*' A v 



0, 



(61) 



which is the standard constraint satisfied by a massive 
Proca field, in which — R$>(A 2 )/2k plays the role of an 
extra contribution to the vector's potential V(A 2 ). 

The Einstein and conservation equations, obtained re- 
spectively by varying with respect to the Jordan metric 
and the matter fields, yield 

^{A 2 )G^ - (V M V„ - g^D) *(A 2 ) + m'{A 2 )A^A u 



F^F v a - -g^F 2 + 2V'A»A V - V 9liV + T™ 



V T^ v — 



(62) 
(63) 



the second equation being no surprise since the matter 
fields are minimally coupled to the Jordan metric. 

On the other hand, Eq. (fB"2"]) already exhibits the 
deadly problem that this class of models presents: Some 
gauge-dependent second derivatives of the vector field A^ 
are generated in the left hand side. They come from the 
Riff (A 2 ) term in action |59|) . which breaks the gauge in- 
variance of the vector's kinetic term. Indeed, the scalar 
curvature R contains second derivatives of the metric, 
therefore, after integration by parts, second derivatives 
of A^ which cannot be written in terms of the gauge- 
invariant Faraday tensor _F),„ (nor its dual F^ v ). We 
thus expect to excite the generic helicity-0 ghost of non- 
gauge-invariant vector theories, as in Sec. Ill Al above. We 
will see below that this will become explicit thanks to a 
change of variables, namely by rewriting the same the- 
ory in the so-called Einstein frame. Equation (|6"2"|) also 
illustrates why this ghost is never noticed when studying 



Indeed, the gauge-dependent second derivatives are act- 
ing on a function of A 2 , and therefore disappear at linear 
order in A^. This is actually already manifest in action 
(|59[) . since the gauge-dependent terms involving deriva- 
tives of A^ are of the cubic form A 2 ddh (where h denotes 
schematically a perturbation of the metric), and there- 
fore of quadratic order in the field equations. 



2. Einstein frame 

The kinetic terms of the spin-1 and spin-2 degrees of 
freedom are not diagonalized in action (|59")) , as clearly il- 
lustrated by the field equations (f6"D")K (|62 p . As for scalar- 
tensor theories, the theory is better analyzed in the so- 
called Einstein frame, defined by diagonalizing the ki- 
netic terms. This can be achieved thanks to a conformal 
rescaling of the metric 

9 ; u = y(A 2 )g,„. (64) 

For the sake of clarity, we set A 2 = g* v A^A V so that 

A 2 



A: 



^>{A 2 y 



(65) 



which is assumed to be invertible as a function ${A 2 ) = 
A 2 . When performing the conformal transformation and 
also replacing A 2 in terms of A 2 , we obtain that ac- 
tion (1551) can be rewritten as 



5' 




R* 



jLz 2 {Al){d»Al) 2 -\f 2 

4/C 4 



^matter [tpml B (A^) g* v ] 



(66) 



where only use of the Einstein metric g* is made in 
all contractions and in defining the Ricci scalar We 
notably define as usual F^ u = d^A v — d u A^ but Ft" = 
g* P 9TF p(7 and At = gt ' A v . The three functions of A 2 
that appear in this action are given by 



B{Al) 
Z{Al) 
W(Al) 



1/*(A 2 ), 
dlnB 



^{A 2 )^'{A 2 ) 



dA 2 ^(A 2 ) 
V{A 2 )/^ 2 {A 2 ). 



A 2 W{A 2 



(67) 
(68) 
(69) 



linear perturbations, around a background where A^ 



0. 



The kinetic terms of the vector A p and the tensor g* v 
are now diagonalized in action (|66p . in a covariant way. 
This will allow us to consider the vector sector alone in 
in Sec. IIII A 31 below, say in a freely falling elevator, to 
analyze its stability. 

Let us however underline a subtlety related to vector 
fields in curved spacetime, as soon as their kinetic term 
is not a mere function of the Faraday tensor F^ v and 
its dual F^ v . Indeed, the contribution proportional to 
Z 2 in action (|6"6"|) involves a cross kinetic term of the 
form dAdg* , because the inverse metric enters the square 
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Al = g'i^AaAp. This can be seen either by writing 
d^A 2 = 2A^d p A a + AaApd^g^ in a non-covariant way, 
or by recalling the presence of a Christoffel symbol in 
the covariant form d u A 2 



This is also il- 



lustrated by the Einstein equations deriving from action 
(luT)l) . which read 



G* = k (T. 



*mat 



//// 



■T, 



•■EM 



2\ 2 



Z^^Ald^Al 



zz> (d a Aiy + z 2 n*Al 



An.Au 



';Z 2 (d a Al 



9 LIU 1 



(70) 



where 



{1/ \/-g*){5S ma .t/ Sg*^) is the matter 
energy-momentum tensor as defined in the Einstein 
frame. The presence of second derivatives of the vec- 
tor field in Eq. ([70]) . in the form of D*A 2 , underlines 
that cross kinetic terms were actually still involved in ac- 
tion (|66[) . On the other hand, no curvature tensor enters 
the Maxwell equations deriving from action ([66)) in the 
Einstein frame: 



v;Fr 



a: 



2W' - ~ZZ' {d a Al 



—z 2 n*Al 



(71) 



where T* 



9p.u mat * 



It should be noted that the 



actual energy-momentum tensor measured by an ob- 
server is the Jordan-frame one, defined as T^ at = 
(2/ ^/~ r g)(5S ma . t /8g pu ), and its trace as T mat = g^T^. 
It is related to its Einstein-frame counterpart in a non- 
trivial way, because B(A 2 )g* tlv depends on the Einstein 

metric g pu also through A 2 = g*A a A@. One finds 
r mat = B3 { T Zt ~ B'A^A^T^t), so that the last term 
within the square brackets of Eq. (|7ip may also be writ- 
tenas* 2 *'T mat = ZT mat . 

Although the kinetic terms are covariantly diagonal- 
ized in the Einstein- frame action (|66[) . one may be wor- 
ried by the non-covariant cross term dAdg* it still con- 
tains. Indeed, it is well known that such cross terms 
may contribute positively to the kinetic energy of a de- 
gree of freedom. The best known example is Brans- 
Dicke scalar-tensor theory, defined by the action S = 
J d 4 xy^g [$R - (w/$)(9 M $) 2 ] , where the spin-0 degree 
of freedom carries positive energy provided lo > — | . For 
— | < uo < 0, one may thus naively think the scalar field is 
a ghost, but the cross kinetic term involved in <I>i? (after 
partial integration) is enough to guarantee the positivity 
of energy. To check that the remaining cross kinetic term 
of action ([66]) actually does not change our conclusion of 
Sec. IIII A 31 below, let us eliminate it in a non-covariant 
way. The clearest way to do so will be to start again 
from the Jordan- frame action (|59[) . and to consider per- 
turbations around a given background, keeping covariant 
expressions with respect to the background metric. Let 
us define g^ 1 = g^ v + and A^ 11 = A^ + a^, and ex- 
pand (|59p to second order in the dynamic perturbations 



h^ v and a M , using the background metric g^ to contract 
indices or define covariant derivatives. The kinetic terms 
of these perturbations then read 

-±-^{A 2 )v,h a(j (2 3 Q V 5 - g ap g 1& ) V/m 

+ T ^(^)(2V^-V,/ i ) 2 

_ -L*'(A 2 )(V„^ - V p h)(2A p V p a p - A p A°V^h pa ) 

ZK 



--(V M a„ - V^) 2 , 



(72) 



where h = g a ^h a p is the trace of the Jordan metric per- 
turbation. The first two terms of (jT2[) are the standard 
kinetic term of a spin-2 graviton, multiplied by a global 
factor iff (A 2 ) depending on the background vector field 
A^, the fourth term is the standard Maxwell kinetic term, 
and the third term exhibits the cross kinetic terms V/iVa 
generated by the nonminimal coupling R^(A 2 ) of action 
(|59"|) . Before diagonalizing these kinetic terms, let us re- 
call that the general coordinate-invariance of action ([59)) 
implies the gauge-invariance of the Jordan metric per- 
turbation hap (although the Jordan metric g a does not 
describe a pure spin-2 degree of freedom). We may thus 
fix the harmonic gauge in Eq. (|T2[) by imposing 

2V v h v ll = V^h. (73) 

This choice not only removes the second term of (|72jl. 
but also simplifies the third term as 



V(A 2 ) 
4k 



Wh (2A p V^a p - AM CT V M V)' ( 74 ) 



It is now straightforward to check that the redefinition 



incw 
rl afd 



2f'A p a n 



(*- A 2 *') 2 



[(* + A 2 V') 



9afi 



-2(A 2 W) 2 
WA a Ap] 



(75) 



then suffices to eliminate all cross terms V/i now Va. This 
change of variable differs in several ways from the confor- 
mal transformation (|64[) used above in the covariant cal- 
culation. Indeed, it now contains a "disformal" (i.e., non- 
conformal) contribution proportional to A a Ap. More- 
over, it clearly breaks general covariance since the mod- 
ification of h a /3 is proportional to the mere contraction 
A p a p , whereas the expansion of A 2 ull — A 2 + 2A p a p — 
A p A"h pa + ((h, a) 2 ) also involves the projection of the 
metric perturbation along the background vector field, 
A p A a h pa . This is one of the reasons why (f75)) allows 
us to cancel the cross kinetic term dAdg* we had in the 
covariant action ([66]) . Finally, Eq. ([65)) happened not 
to be invertible for the simple case of ty(A 2 ) — A 2 , as 
also illustrated by the vanishing denominator of defini- 
tion (p5)) . whereas Eq. (|75p is always invertible as soon 
as V(A 2 ) f 0. 

It should be noted that the gauge fixing ([73")) is non 
trivial in terms of the new variable /i^ w , since it now 
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also involves the vector perturbation a p . However, as 
already underlined above, the general covariance of the 
Jordan-frame action ([59]) anyway guarantees this choice 
is allowed. It underlines that some cross-kinetic terms 
(between h ncw and a) are actually pure gauge, and cannot 
contribute to any physical observable. Replacing now 
definition (|75|) in (f72")) . still in the gauge (|75|) . we can 
read off the full kinetic term of the vector perturbation: 



Note that at linear order in the field equations (i.e., 
quadratic order in the action or the Hamiltonian), we re- 
cover 7T° = as in gauge-invariant vector theories. There- 
fore the ghost instability present in the nonminimally 
coupled models ([59]) or (|66[) cannot be noticed when 
studying first-order perturbations (around a vanishing- 
vector background). 

We deduce that the Hamiltonian density takes the form 



1 



(V^cv-V^) 2 - 



k (* - A 2 ^ 1 ) 2 + 2{A 2 ^') 



(76) 



This is similar to the expression (J66J) we found in the fully 
covariant case, with the minor difference of a global fac- 
tor for the second term [coming from the fact that 
we use the Einstein metric (|64l) to contract all indices in 
(|6"6")) . whereas we kept the original Jordan metric g pv as 
our present background], the important difference that 
the denominator of this second term contains a contri- 
bution +2(A 2 ^') 2 in addition to the square - A 2 ^') 2 
coming from Z 2 [this change also comes with a modifica- 
tion of the global numerical factor from 3 to 2], and the 
crucial difference that all cross kinetic terms have been 
cancelled. When considering (|6"6"|) in a flat background 
g* = ri^v (or in a Fermi coordinate system), we thus 
get an expression of the same form as expansion (|76p . 
the only difference being the precise definition of Z . In 
Sec. IIII A 31 below, the nonvanishing of this function Z 
will be the only needed information, therefore one may 
work with the covariant action (|66[) although its kinetic 
terms are not fully diagonalized. 



Both (|66|) and ([76)) show that the mode of a p which is 
polarized in the direction of the background A p behaves 
as if it were a positive-energy scalar field (see also Sec. 5 of 
the recent Ref. [Hj]). However, it is coupled to the other 
vectorial modes via the standard Maxwell kinetic term, 
and we will see now that this causes a deadly instability 
of the model. 



3. Hamiltonian analysis 

The stability analysis of any model is much more easily 
performed in the Einstein frame, where the spin 2 and 
the other degrees of freedom decouple. As discussed in 
the previous section, there still exists a cross kinetic term 
dAdg* in the covariant action (|66p . but eliminating it in a 
non-covariant way, as in Eq. (|T6[) , keeps the same general 
form for the vector's kinetic term. Let us thus consider an 
action of the form (|6"6"|) , with Z ^ but maybe different 
from (f6"8")h and focus on the vector's dynamics in a flat 
geometry g pi/ = r\ pv . The conjugate momenta are then 
given by 



--Z 2 A Q d t (Al), (77) 

K 

A l ~d l A Q + -Z 2 A l d t (Al). (78) 



H = 



l -F 2 

4 v 



W 



— Z 2 x \d l (Al)] 2 + — ( 

4k 1 V * ;J 12 \ZA 



-0 



\ (< + T Ai + diAo 



\ (drAof 



(79) 



Since ir° is not identically zero in Eq. (|77[) . the Aq com- 
ponent is dynamical and independent from the spatial 
components Ai . We may thus consider a particular back- 
ground such that Aq ^ while 7r° = 0, Ai = and 
it 1 = —diAo, and the Hamiltonian density then reads 



n ~ w 



(cUo) 



(80) 



Initial data of the form Aq — esin(a:/e 2 ), with e — > 0, 
would thus make this Hamiltonian density tend towards 
— oo. This suffices to show that the nonminimally coupled 
vector model (15^1) or (16U1) is unstable. 



F case 



B. 



We can also consider theories in which the Faraday 
tensor is nonminimally coupled to the Ricci scalar. In 
the Jordan frame, such theories will have an action of 
the form 



S = J d^V^g 

H~ ^matter [V^m ; Q\xv 



^S(F 2 ) - f(F 2 ) V(A 2 ) 



(81) 



with the same definitions as in the previous sections. 

One may be tempted to introduce the analogue of an 
Einstein metric by defining 



9. 



S(F 2 )^, 



(82) 



but since this definition involves derivatives of the vector 
field, it cannot be used consistently in a Lagrangian (see 
footnote [3] above) . 

Actually, because the scalar curvature R involves sec- 
ond derivatives of the metric tensor g^ „, action (|8ip gen- 
erates third derivatives of the vector field in the metric 
field equation, and third derivatives of the metric (i.e., co- 
variant derivatives of the curvature tensor) in the vector 
field equation. Initial data on a Cauchy surface should 
thus contain more information than the values of the 
fields and their time derivatives. Therefore, this class 
of models must involve some extra degrees of freedom, 
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in addition to the vector and the metric we wished to 
introduce. Such higher derivatives are known to pro- 
duce generically ghost degrees of freedom, i.e., to cause 
the theory to be unstable. This is notably the conse- 
quence of a theorem by Ostrogradski [48[ , well discussed 
in Ref. [39j ] . However, this theorem can be applied only 
on so-called "nondegenerate" Lagrangians, which pro- 
duce fourth-order field equations. Therefore, we are here 
in a typical case where we expect a serious instability to 
manifest, but where we cannot use the generic theorem 
which proves so without any ambiguity. 

To understand intuitively the instability of a theory 
defined by action (|81[) . one may consider a toy model 
involving two coupled scalar fields in flat spacetime, 
C = -(d^ip) 2 - {d^ip) 2 + \{d^) 2 {d v il}) 2 , where A is a 
coupling constant (see Sec. V A of Ref. [13 )• Here ip 
and i/j play the roles of the metric tensor and of the vec- 
tor field of Eq. (|81[) . The corresponding Hamiltonian 
density reads H = p 2 + ijj 2 + {d l p) 2 + (d^) 2 + 4\tp 2 ip 2 - 
\[p 2 + (dip) 2 ]{ip + (djtp) 2 ], and it can be made arbi- 
trary large and negative whatever the sign of A. Indeed, 
if A < 0, it suffices to choose a homogeneous configu- 
ration di p = di%jj = and large enough values of ip 2 
and ip 2 . On the other hand, if A > 0, instantaneously 
constant fields tp — tp — with large enough spatial 
derivatives (dip) 2 and (diip) 2 suffice to make TC tend 
towards — oo. Even more intuitively, in a given back- 
ground of ip, the second scalar field ip behaves as if its 
kinetic term were multiplied by [1 — X(d^p) 2 ]. If the 
^-background is chosen such that X(d tl p) 2 be negative 
enough, then ip will behave as a ghost, and its contri- 
bution to the Hamiltonian density will be unbounded by 
below. Therefore, there do exist field configurations such 
that 7i is as negative as one wishes, and this proves the 
instability of the toy model. Such a hand-waving argu- 
ment can now also be used on more involved models, 
for instance C = — {d^p) 2 — {d^ip) 2 4- XDp(d t _ l ip) 2 , which 
looks a little more like Eq. (f8"Tj) . where Dp plays the role 
of the scalar curvature R, involving second derivatives. 
The Hamiltonian analysis is now much more involved, 
because the presence of third derivatives of the fields in 
their equations implies the existence of new excitations 
(and the standard Ostrogradski definition of conjugate 
momenta cannot be followed because we are in a degen- 
erate case). But it is still clear that in a given back- 
ground where \Dp is large enough, then ip behaves as 
a ghost and can make the Hamiltonian density tend to- 
wards — oo. Now, if we try to apply this argument to 
action (fSTj) itself, we understand that we need to con- 
sider large enough (positive or negative) curvatures R 
such that the nonminimal coupling RF 2 could change 
the global sign of the vector kinetic term. Particular 
case might thus be safe, for instance if one needs to be 
in the interior of a black hole horizon to reach such a 
condition. Moreover, one may devise models such that 
the function S(x) = £o + Ox + £,2X 2 + ■ • ■ does not con- 
tain any linear term. Therefore, the above hand-waving 
argument does not prove that all models (|81|) are unsta- 



ble, although we do expect so because of the presence of 
higher derivatives in their field equations. We will any- 
way disregard this class of models, because such higher 
derivatives mean that they involve extra degrees of free- 
dom, in addition to the single spin-1 and spin-2 fields we 
wished to consider. 



IV. DIMENSIONAL REDUCTION OF 

LOVELOCK INVARIANTS AND 
COSMOLOGICAL PHENOMENOLOGY 

A Lovelock invariant is defined in even dimen- 
sion D as a Lagrangian density proportional to = 

involving thus a product of D/2 Riemann curvature 
tensors. The best known examples are the cosmological 
constant A corresponding to D = 0, the Einstein- Hilbcrt 
Lagrangian R corresponding to D = 2, and the Gauss- 
Bonnet density RJt Vpa ~ ^R^ + R 2 corresponding to 
D = 4. The integral of Lr> over a D-dimensional space- 
time gives a number depending only on the topology, 
therefore its variational derivative vanishes and it does 
not contribute to the field equations. In dimensions 
lower than D, the density Lr> vanishes identically. On 
the other hand, defines a nontrivial dynamics when 
considered in dimensions higher than D (like R or A in 
4 dimensions). But in spite of the presence of several 
Riemann tensors (for L04), each of them involving 
second derivatives of the metric, the corresponding 
field equations remain of second order. Indeed, any 
third (or higher) derivative must appear in a form 
similar to -R/ii^^i^w multiplied by the antisymmetric 
Levi-Civita tensor gW* 3 — and therefore vanishes by 
virtue of the Bianchi identity Rv^^i^;^] = 0. The 
absence of higher-order derivatives in the field equations 
does not guarantee the stability of the corresponding 
models, but it proves at least that no extra degree of 
freedom is excited, and that the generic ghost modes of 
higher-order theories are avoided. If the Gauss-Bonnet 
density i?^ l/p(T — 4i?^ y -I- R 2 is considered in 5 dimensions, 
for instance, it does contribute to the field equations, 
but keeping them of second order. When performing 
a Kaluza-Klein dimensional reduction, where is 
interpreted as a vector field in four dimensions, we 
thus get a nontrivial vector-curvature coupling which 
does not generate higher-order field equations, and 
avoids thus the deadly instabilities caused by ghost 
modes 7 . We will analyze below the cosmology generated 
by such a coupling. Similar models can be constructed 
by considering the dimensional reduction of higher-order 
Lovelock invariants L§,L%, . . . , and even more general 



Moreover, the dimensional reduction of Lovelock invariants al- 
ways generates (gauge-invariant) combinations of the Faraday 
tensor F^, therefore the ghostlike mode Ao is never excited ei- 
ther; see Sec. Ill Al 
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vector models coupled to both curvature and scalar fields 
are obtained by dimensionally reducing the so-called 
"Galileon" models recently introduced in Ref. [HI and 
generalized in curved spacetimes in Refs. (50[. As we 
will see below, even the simplest case of a dimensionally- 
reduced Gauss-Bonnet density L4 suffices to generate an 
interesting cosmological evolution for the vector field. 



A. Nonminimal couplings to the Riemann tensor 



B. Cosmological dynamics 



We consider the class of models 



5' 



+ ^matter [4^m j 9llv] 



(83) 



with the same notation as in the previous sections, and 
where £, 77 and C are constant parameters. Such theories 
lead to generalization of the Maxwell theory that imply 
a variable speed of light [5l|, [HI, [H, [54| (i-e. propaga- 
tion velocity of the vector field if identified to the one 
describing the photon [EH). However, as shown in [H3 |. 
the corresponding field equations are of second order if 
and only if the parameters £, 77 and £ satisfy 



77 + 2£ = 0, C = £, 



(84) 



and this is precisely what is obtained by dimensionally 
reducing the Gauss-Bonnet density L4 written in a 5- 
dimensional spacetime [57], [H| ■ This can also be checked 
explicitly by deriving the vector field equations 

(1 - £,R)F<". V - 77 (i?^ A % - R V X F X %) 
1 



<R w Fpa . v _ 
(77 + 20^.^ = 0, 



~(2Z + r})R u F^ 



(85) 



in which third derivatives of the metric occur (in the 
form of first derivatives of the curvature tensor) unless 
relations ([84")) are satisfied. Similarly, the Einstein equa- 



tions involve third derivatives of the vector A„ unless 



are satisfied. Since such higher derivatives would 
excite new, generically ghostlike, degrees of freedom, im- 
plying the instability of the model, we will restrict our 
study to the particular case ([84")) . However, second-order 
field equations do not suffice to warrant the consistency 
of the model. These equations should also be hyperbolic, 
and the corresponding Hamiltonian should be bounded 
by below. We will not perform here this analysis, be- 
cause it is even more complex than in the case of cou- 
plings like R^A^A". However, we wish to emphasize 
that this particular class of models offers an interesting 
phenomenology for cosmology and should thus deserve 
more attention. 



Consider action ([83)) where the matter fields reduce 
e.g. to a single scalar field <fi evolving in a potential 
v(<p) that is assumed to drive an inflationary phase in 
the early universe. We consider the vector field as a test 
field whose evolution is then given by Eq. ([83)). Using the 
same notation ([23)) as in Sec. Ill B 31 above, homogeneity 
in a Friedmann-Lemaitre spacetime with metric ([1} im- 
plies diA^ — 0, so that the only nonvanishing component 
of the Faraday tensor is 



F Ql = a{Bi + HB % ) 



(86) 



We recall that the Weyl tensor of a Friedmann-Lemaitre 
spacetime strictly vanishes, C^ vpa = 0, and that (re- 
stricting to a spatially Euclidean spacetime) the nonva- 
nishing components of the Ricci tensor are given by 

R%=3(H + H 2 ), R) = (H + 3H 2 )6], (87) 

so that the only nonvanishing components of the Rie- 
mann tensor are (see e.g. Ref. [59]) 

R l ]ml = a 2 H 2 (S l ml3l - 8\ ljm ), R° i0 j = a 2 (H + H 2 ) lt] . 



The evolution equation for Bi then reduces to 

1 - 6£(H + 2H 2 ) - 7 ? (4iJ + QH 2 ) - 2((H + H 2 ) 
x (F l ° + 3HF ia ) + -(6£ + 47/ + 2()(H + AHH) 
+4(?7 + ()HH F z0 = 0. (89) 

Restricting to the conditions (I84|) , it leads to the equation 

'2 ■ 



(1 + ^)^ + 3 



l +ri (-H + H< 



HBi 



[l + 3rjH 2 ) H + 2(1 + ?]H 2 ) H 2 B, = 0. (90) 



Let us first assume that the universe is undergoing a 
slow-roll inflationary phase close to a de Sitter phase, so 
that we can assume H ~ const, and —H/H 2 = e <C 1 
(e > in most slow-roll inflationary models). If the pa- 
rameter 77 is chosen to be negative, then a fine-tuned 
value H 2 w —(l + s)/r] is such that Eq. ([9"0)l reads 
Bi + (1 - 2e)HB l - 3eH 2 B l = 0, and therefore does 
not involve any undifferentiated Bi at lowest order in e 
[this can easily be made exact thanks to an even finer 
tuning of H(i)]. The two solutions of this equation are 
thus a decaying mode Bi ~ exp[— Ht] ~ 1/a and an al- 
most constant one Bi ~ exp[3eHt] ~ a 3e — even slightly 
increasing if e > 0. It follows that a slow-rolling vector 
field can survive the expansion, contrary to the standard 
lore on vector fields, but at the price of a fine tuning of 
the expansion rate H, related to the nonminimal vector- 
gravity coupling constant 77. 

Let us also consider the dynamics of the vector field 
assuming the background dynamics is given by a(t) cx t p 
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(p = 1/2 for a radiation-dominated universe, p — 2/3 
for a matter-dominated universe, and the limit p » 1 
corresponds to a power-law inflationary model with e = 
1/p). Equation ([50)) then reduces to 



(1 + V H 2 )B\ 



2-i 



3 - - )nH z 

V 



HBi 



2--UE- 



H Bi = 0. 



(91) 



In the case of inflation, we find again that the field is 
diluted unless one imposes the previous fine tuning 1 + 
i]H 2 = — 1/p, which requires r\ < 0. 

To discuss the dynamics during the matter and 
radiation-dominated era, let us introduce the time scale 
r * = P\fV\\- in t ne radiation era, the coefficient of Bi 
is always proportional to r\ (instead of being zero in the 
standard case). At early times (t <C r*), Eq. ([91) re- 
duces to Bi — Bi/{2t) — Bi/t 2 — which has two so- 
lutions, a decaying mode oc \j\ft oc a~ 4 and a grow- 
ing mode oc t 2 oc a 4 while, Ci = Bi + HBi behaves as 
d oc t oc a 2 . At later times (t » r*), Eq. ([9Tj) reduces 
to B, + 3Bi/(2t) - r/_B,;/(4i 4 ) = 0, which differs from the 
standard equation by the term proportional to Bi. The 
solutions of such an equation are given in terms of Bessel 
functions and will be oscillating if rj < while they have 
a mode oc (t / r*) -1 / 4 i^/^-r* /t) if r\ > that grows and 
then freezes to a constant. These behaviors at early and 
late times differ from the standard dynamics of a vec- 
tor field and exist whatever the value of rj. In the matter 
era, the dynamics is only modified at early times (t <C 7%) 
since Eq. reduces to B,+9tB l /(2T])-10B t /(9t 2 ) = 0, 
the main modification arising from the fact that the co- 
efficient of Bi is now proportional to l/H and not to H 
anymore. This equation has a growing mode. 



C. Discussion 

In this class of theories, a slow-rolling vector field can 
survive during inflation, contrary to the standard lore on 
vector fields, but at the price of a fine tuning of the ex- 
pansion rate H, related to the nonminimal vector-gravity 
coupling constant 77. It requires that 77 be negative and is 
related to the energy scale of inflation by \r\\ ~ 1/Hf n ^. 
The general action should thus contain terms of the form 



1 



1 



1 



C D -M;R - -F - -r]H z F ( 



0( • 



where M p is the Planck mass. Then, since during infla- 
tion R ~ l2H 2 n j while we need i]H 2 n j ~ — 1, the correc- 
tion term to the standard Einstein-Maxwell Lagrangian is 
of the order of F 2 J2 = -F 2 /A. As long as H mf /M p < 1, 
as is usually the case in inflation, this correction is negli- 
gible compared to the Einstein-Hilbert while being of the 
same order as the Maxwell term. In spirit, this solution 
leading to a slow-rolling vector field is similar to the one 



invoked in Refs. [1,[T3|, which used a coupling of the form 
£RA 2 , that we saw to be unstable. We cannot prove at 
this stage that this will not be the fate of this model that 
needs to be analyzed in detail. 

We have also seen that the dynamics during the radia- 
tion and matter-dominated eras allows for growing solu- 
tions whatever the value of the parameter 77. This opens 
an interesting phenomenology that we postpone to fur- 
ther study. 



V. CONCLUSIONS 

In this article, we have investigated general models of 
vector fields that have recently been considered in cos- 
mology, in relation with a source of anisotropy or the 
construction of MOND-inspired field theories. 

We have shown that the class of /(F 2 )-theories suf- 
fers from hyperbolicity problems, while both f(F 2 ) and 
f(F 2 ,FF) models predict a dilution of the vector field 
during the cosmological expansion. 

When allowing for a nonminimal coupling to the met- 
ric, we have proven that the class of /(^4 2 )i?-theories has 
a Hamiltonian which is unbounded from below, while the 
f(F 2 )R- models involve higher derivatives of the fields 
and thus contain extra degrees of freedom (which are 
generically expected to carry negative energy). 

These results set strong constraints on vector field 
models, as long as they are considered as fundamental 
theories — i.e., notably, that no field entering the ac- 
tion is considered as a fixed background that cannot be 
varied. [From a theoretical point of view, let us remind 
that an action is not just a list of symbols but involves 
also the definitions of these symbols, e.g. what are the 
fundamental fields; see the discussion of the difference be- 
tween and d^cj), or the difference between a potential 
and a Lagrange multiplier.] But even as effective models, 
the constraints we derived for their stability and causal- 
ity should always be satisfied in their domain of validity, 
and at least in the domain where their cosmological evo- 
lution is studied. It happens that to avoid the dilution of 
the vector field during the expansion of the Universe, one 
would need the nonlinear terms to be of the same order of 
magnitude as the main kinetic term, i.e., precisely in con- 
ditions where the positivity of the Hamiltonian and the 
well-posedness of the Cauchy problem should be checked 
carefully. 

To finish, we pointed out that in the class of theo- 
ries obtained by dimensional reduction of Lovelock in- 
variants, there exist cases that allow for the existence of 
a slow-rolling vector field. Although we did not study 
the boundedness by below of the Hamiltonian nor the 
hyperbolicity of the field equations, because of their com- 
plexity, we underlined that the field equations remain of 
second order in spite of the nonminimal coupling of the 
vector field to curvature. Such models contain thus only 
the spin-1 and spin- 2 degrees of freedom we wished to 
consider (in addition to other matter fields), and they 
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are phenomenologically quite appealing for cosmology. 

Similar models as the one we studied in Sec. lIVI are ob- 
tained by dimensionally reducing higher-order Lovelock 
invariants, and more general tensor-vector-scalar models 
yielding second-order field equations can also be defined 
by dimensional reduction of Galileon actions [i^,[53| writ- 
ten in more than 4 dimensions. It is also possible that 
the construction of scalar Galileons can be generalized to 
vector fields, yielding nonminimal vector-curvature cou- 
plings of a different nature than those obtained from 



Lovelock invariants. All such second-order models de- 
serve being studied both mathematically and for their 
phenomenological predictions in a cosmological context. 
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